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Abstract 

Preliminary  results  for  numerical  simulations  of  regional  propagation  in  realistic  earth  models 
are  presented.  The  modeling  is  motivated  in  part  by  three-component  array  observations  which 
show  that  the  P  and  5  codas,  and  Lg  wavetrains  for  regional  sources  axe  comprised  of  substantial 
amounts  of  forward  scattered/multipathed  energy,  and  in  part  by  results  from  numerical  simula¬ 
tions  of  teleseismic  wave  propagation  which  suggest  a  correlation  between  the  coda  level  and  rate 
of  decay  and  the  aspect  ratio  of  lithospheric  heterogeneities.  Numerical  simulations  are  conducted 
using  a  two-dimensional,  fourth  order  accurate,  P-SV  finite  difference  routine.  The  models  con¬ 
sist  of  a  heterogeneous  layer  over  a  homogeneous  half-space.  Perturbations  are  in  velocity  and 
parameterized  using  a  gaussian  correlation  function. 

The  wavefields  generated  using  models  containing  high  aspect  ratio  heterogeneities  (AN  ISO)  ex¬ 
hibit  several  characteristic  consistent  with  regional  array  observations.  First,  the  Pg  wavetrain  and 
coda  consists  of  forward  scattered,  coherent  plane  wave  energy  trapped  in  the  crustal  waveguide. 
Models  containing  isotropic  heterogeneities  (ISO)  inhibit  the  extent  to  which  energy  is  trapped  in 
the  crustal  waveguide  and  produce  codas  dominated  by  randomly  scattered  energy.  Waveforms 
generated  using  models  containing  ANISO  heterogeneities  exhibit  another  feature  commonly  ob¬ 
served  in  regional  array  data.  That  is,  that  detailed  features  in  waveforms  observed  at  adjacent 
sensors  exhibit  a  high  degree  of  coherence,  while  those  observed  at  sensors  separated  by  distances 
on  the  order  of  a  wavelength  exhibit  a  somewhat  surprising  loss  of  coherence.  Waveforms  gener¬ 
ated  using  models  containing  ISO  heterogeneities  exhibit  a  degree  of  incoherence  inconsistent  with 
regional  array  observations. 

The  modeling  results  presented  in  this  report  add  to  the  growing  body  of  evidence  illustrating 
the  important  role  spatially  anisotropic  heterogeneities  play  in  seismic  wave  propagation  and  scat¬ 
tering.  An  important  but  overlooked  implication  inherent  in  parameterizing  heterogeneities  with 
spatially  anisotropic  correlation  functions  is  that  the  dimensionless  normalized  wavenmnber  ka 
commonly  used  to  categorize  the  scattering  regime  looses  its  meaning.  Estimates  of  the  effective 
scattering  cross  section  and  the  associated  scattering  radiation  pattern  need  include  the  vector 
quantity  k-a  to  reflect  the  fact  that  the  nature  of  scattering  depends  not  only  on  the  (scalar)  ratio 
of  the  "size”  of  the  heterogeneities  to  the  wavelength  of  the  incident  wavefield,  but  also  on  the 
spatial  relation/orientation  of  the  incident  wavefield  and  heterogeneities. 
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Objective 


The  main  objective  of  our  research  is  to  provide  an  improved  understanding  of  regional  wave 
propagation  and  to  aid  in  the  characterization  of  lithospheric  heterogeneities  by  combining  analysis 
of  observed  three-component  array  data  and  numerical  simulations  of  elastic  wave  propagation  in 
realistic  earth  structures. 


•  Preliminary  Research  Results 


This  report  provides  prehminary  results  from  numerical  simulations  of  regional  wave  propagation  in 
realistic  earth  structures.  The  modehng  is  motivated  in  part  by  three-component  array  observations 
which  show  that  the  P  and  S  codas,  and  Lg  wavetrains  for  regional  sources  are  comprised  of  a 
considerable  amount  of  forward  scattered/multipathed  energy  (Wagner  &  Owens,  1993,1995a,b; 
Dainty  &  Schultz,  1995),  and  in  part  by  results  from  numerical  simulations  of  teleseismic  wave 
propagation  which  suggest  a  correlation  between  the  coda  level  and  rate  of  decay  and  the  aspect 
ratio  of  lithospheric  heterogeneities  (Wagner  &  Langston,  1992).  The  results  presented  in  this 
report  are  simply  an  extension  of  the  teleseismic  modeling  to  the  regional  problem  in  an  effort  to 
determine  the  types  of  structure  required  to  produce  wavefields  that  exhibit  the  same  characteristics 
observed  in  regional  array  data. 

The  basic  model  used  in  our  investigation  is  a  heterogeneous  layer  (i.e.  Mthosphere)  over  a 
homogeneous  half-space  (i.e.  mantle).  The  use  of  this  general  model  for  both  teleseismic  and 
regional  simulations  is  motivated  by  the  assumption  that  small  scale  heterogeneities  are  generally 
confined  to  the  hthosphere.  The  heterogeneities  are  parameterized  using  a  gaussian  correlation 
function  (space  domain  wavenumber  domain) 
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where  Oj,  are  the  horizontal  and  vertical  correlation  distances,  respectively.  The  two  end 
members  of  this  model  parameterization  scheme  are  models  with  isotropic  heterogeneities  (ISO)  and 
models  with  homogeneous  plane  layers.  As  the  aspect  ratio  (oj/oz)  increases,  the  heterogeneities 
are  stretched  in  the  horizontal  direction.  Heterogeneities  with  higher  aspect  ratios  (ANISO)  mimic 
plane  layers  of  varying  lateral  extent.  These  types  of  models  are  intuitively  appealing  because 
they  agree  with  geologic  models  where  layers  do  not  have  infinite  lateral  extent  but  are  truncated 
by  faulting,  pinch  outs,  etc..  Figiure  1  shows  examples  of  ISO  and  ANISO  models;  the  correlation 
distances  are  ax=az=5  km  for  the  ISO  model,  and  Ox=50  km,  az=5  km  for  the  ANISO  model. 

The  results  presented  in  this  report  are  for  models  whose  heterogeneities  consist  of  correlated 
P  and  S  velocity  pertmbations.  The  mean  P  and  S  velocities  are,  respectively,  5.25  and  3.03  km/s 
in  the  layer,  and  7.50  and  4.33  km/s  in  the  half-space.  The  standard  deviation  of  the  P  and  S 
velocity  perturbations  is  5%  of  their  respective  backgroimd  means.  The  density  is  2.5  kg/m^  for 
both  the  layer  and  the  half-space.  No  intrinsic  attenuation  is  included  in  the  model. 

Numerical  methods  are  currently  the  only  methods  available  to  compute  waveforms  in  a  medium 
whose  properties  vary  both  laterally  and  verticedly  in  a  random  manner.  Simulations  were  per¬ 
formed  using  a  two-dimensional,  fourth-order  accurate,  displacement,  P-SV  finite  difference  (FD) 
routine.  Because  our  main  objective  was  to  identify  models  that  will  produce  predominantly  for¬ 
ward  as  opposed  to  back  scattered  codas,  it  was  crucial  that  we  use  a  numerical  approach  capable 
of  modehng  both  the  forward  and  back-scattered  components  of  the  wavefield.  Absorbing  bound¬ 
ary  conditions  (Clayton  &  Engquist,  1977)  are  used  on  the  sides  and  bottom  of  the  model,  and 
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Figure  1:  Examples  of  models  with  spatially  isotropic  (ISO)  and  spatially  anisotropic  (ANISO)  het¬ 
erogeneities.  The  perturbations  are  parameterized  using  a  gaussian  correlation  function.  The  ISO 
model  (top)  has  spatially  isotropic  heterogeneities  with  a  correlation  distance  of  5  km.  The  ANISO 
model  (bottom)  has  spatially  anisotropic  heterogeneities  with  vertical  and  horizontal  correlation 
distances  of  5  km  and  50  km,  respectively. 

stress-free  boundary  conditions  (Ilan  &  Loewenthal,  1976)  at  the  surface.  The  source  is  an  isotropic 
gaussian  at  a  depth  of  5  km. 

Figures  2  and  3  show  the  vertical  and  horizontal  component  snap-shots  for  the  ANISO  and 
ISO  models,  respectively  (note  that  each  frame  is  normalized  to  its  respective  maximum;  time  in 
seconds  is  listed  to  the  right).  While  Figures  2  and  3  do  not  show  propagation  to  regional  distances 
(only  to  125  km  here),  it  is  clear  that  the  salient  wavefield  features  are  established  at  this  early 
stage.  Figures  2  and  3  illustrate  two  important  differences  in  the  wavefields  associated  with  these 
two  models.  First,  the  wavefield  generated  using  the  ANISO  model  is  consistent  with  regional 
array  observations  in  that  it  consists  largely  of  forward  scattered,  coherent  plane  wave  energy 
(note  also  the  coherent  plane  wave  nature  of  the  energy  entering  the  ’"mantle”).  The  wavefield 
generated  using  the  ISO  model  is,  conversely,  dominated  by  randomly  scattered  energy.  The  few- 
plane  waves  generated  in  the  ISO  model  are  the  direct,  and  the  P  and  P-to-SV  converted  waves 
reflected/refracted  at  the  layer /half- space  interface  and  the  free  surface.  While  we  certainly  would 
not  suggest  that  the  ISO  model  represents  a  reasonable  model  for  lithospheric  structure,  it  is  clear 
that  models  with  ISO  heterogeneities  inhibit  the  development  of  wavefields  with  characteristics 
consistent  with  observed  array  data.  The  heterogeneity  aspect  ratio  appears  to  play  the  key  role 
in  determining  the  extent  to  which  energy  is  forward  scattered  and,  consequently,  trapped  in  the 
crustal  waveguide. 
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Figure  2:  Vertical  (left)  and  horizontal  (right)  component  snap-shots  for  propagation  through  the  model  with  spatially  anisotropic  heterogeneities 
(see  Figure  1).  Each  frame  is  normalized  to  its  respective  maximum.  Time,  in  seconds,  is  listed  to  the  right.  Note  that  the  P-coda  is  dominated 
by  coherent,  forward  scattered  plane  waves,  a  characteristic  feature  observed  in  both  local  and  regional  array  data. 
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Figure  3:  Vertical  (left)  and  horizontal  (right)  component  snap-shots  for  propagation  through  the  model  with  spatially  isotropic  heterogeneities 
(see  Figure  1).  Each  frame  is  normalized  to  its  respective  maximum.  Time,  in  seconds,  is  listed  to  the  right.  Note  that  the  P-coda  is  dominated 
by  randomly  scattered  energy,  this  feature  is  inconsistent  with  array  observations  for  both  local  and  regional  sources. 
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Figure  4:  Vertical  component  (velocity)  P„,  Pg  and  P^-coda  waveforms  generated  using  the  ANISO 
(left)  and  ISO  models  (right)  shown  in  Figure  1.  Sensor  spacing  is  1  km.  Source  distance  increases 
from  top  (250  km)  to  bottom  (260  km). 

The  second  significant  difference  illustrated  in  Figures  2  and  3  is  that  the  ANISO  model  facil¬ 
itates  the  excitation  of  considerable  source  zone  reverberations.  This  feature  coupled  with  other 
structural  effects,  such  as  scattering  from  topographic  features,  could  conceivably  result  in  enhanced 
contributions  to  the  regional  wavefield  from  source  zone  scattering  (Greenfield,  1971). 

Figure  4  shows  the  vertical  component  (velocity)  waveforms  generated  using  the  ANISO  (left) 
and  ISO  models  (right)  observed  at  250-260  km.  Because  we  used  a  P  source  we  limit  attention 
to  the  Pn  and  Pg  wavetrains  and  codas.  The  P„  for  both  models  are  similar  and  rather  nonde¬ 
script.  Pn  coda  is  not  readily  produced  in  either  of  the  models,  both  of  which  have  homogeneous 
’’mantles”.  The  implication  here  is  that  coda  is  not  the  product  of  scattering  in  the  crust, 
but  Pn  scattering  that  occurs  in  the  upper  mantle.  Presumably,  therefore,  the  addition  of  ANISO 
heterogeneities  to  the  upper-mantle  would  be  needed  to  generate  forward  scattered  P„  coda. 

The  Pg  wavetrain  and  coda  generated  using  the  ANISO  model  exhibit  several  characteristics 
observed  in  regional  array  data.  First  is  simply  the  extended  coherent  nature  of  Pg.  There  is  no 
one  distinct  Pg  arrival  but  several  equal  amplitude  arrivals,  some  of  which  may  be  a  direct  result  of 
the  aforementioned  source  zone  reverberations,  that  have  traveled  slightly  different  crustal  paths 
from  source  to  receiver  (Kennett,  1989).  The  second  similarity  between  the  ANISO  and  observed 
array  data  is  that  subtle  feature  in  waveforms  observed  at  adjacent  sensors  exhibit  a  high  degi'ee 
of  coherence,  while  those  observed  at  sensors  separated  by  distances  on  the  order  of  a  wavelength 
exhibit  an  unexpected  reduction  in  correlation,  and  those  observed  at  sensors  separated  by  several 
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wavelengths  are  seemingly  unrelated.  The  waveforms  generated  using  the  ISO  model  exhibit  a 
degree  of  sensor-wise  incoherence  inconsistent  with  regional  array  observations.  In  fact,  the  Pg 
wavetrain  generated  using  the  ISO  model  exhibits  almost  no  coherent  coda.  Both  the  extended 
nature  and  slow  spatial  evolution  of  the  AN  ISO  model’s  Pg  wavetrain  result  from  the  interaction 
and  interference  of  the  numerous  forward  scattered  plane  waves  generated  by  scattering  from  the 
spatially  anisotropic  heterogeneities  and  trapped  in  the  crustal  waveguide  and,  potentially,  the 
numerous  smaller  waveguides  bound  by  the  spatially  anisotropic  heterogeneities. 

•  Conclusions  &  Future  Plans 

The  FD  modeling  outlined  in  this  report  w^as  conducted  in  an  effort  to  determine  the  type  of  struc¬ 
tures  capable  of  producing  regional  wavefields  with  characteristics  consistent  with  those  observed 
in  regional  three-component  array  data.  Namely,  wavefields  whose  P  and  S  codas,  and  Lg  wave- 
trains  contain  significant  amounts  of  forward  scattered/multipathed  energy.  The  use  of  models 
with  spatially  anisotropic  heterogeneities  was  motivated  by  previous  results  from  simulations  of 
teleseismic  wave  propagation.  For  the  teleseismic  case,  the  heterogeneity  aspect  ratio  appears  to 
play  an  important  role  in  determining  the  coda  level  and  rate  of  decay  by  controlling  the  extent 
to  which  vertically  propagating  energy  is  scattered  into  the  horizontal  direction.  The  modeling 
results  presented  in  this  report  suggest  that  for  the  regional  case,  the  heterogeneity  aspect  ratio 
plays  an  important  role  in  determining  the  extent  to  which  energy  is  forward  scattered  and,  conse¬ 
quently,  trapped  in  the  crustal  waveguide.  It  is  this  feature,  the  presence  of  waves  such  as  Pg  and 
Lg  comprised  of  energy  trapped  in  the  crustal  waveguide,  that  differentiates  the  regional  wavefield 
from  local  and  teleseismic  wavefields.  Models  with  spatially  anisotropic  heterogeneities  not  only 
produce  wavefields  consistent  with  array  observations,  but  also  provide  a  physically  and  intuitively 
appealing  way  to  parameterize  earth  structure. 

The  modeling  results  presented  in  this  report  add  to  the  growing  body  of  evidence  illustrat¬ 
ing  the  important  role  spatially  anisotropic  heterogeneities  play  in  seismic  wave  propagation  and 
scattering  (cf.  Wagner  &  Langston,  1992;  Hestholm  et  al.,  1994;  Hoshiba,  1995).  An  important 
but  overlooked  implication  inherent  in  parameterizing  heterogeneities  with  spatially  anisotropic 
correlation  functions  is  that  the  dimensionless  normalized  wavenumber  ka  commonly  used  to  cat¬ 
egorize  the  scattering  regime  (i.e.  Rayleigh,  Mie,  Fresnel)  looses  its  meaning.  Clearly,  spatially 
anisotropic  heterogeneities  can  not  be  treated  as  point  scatterers.  The  effective  scattering  cross 
section  and  associated  scattering  radiation  pattern  for  AN  ISO  heterogeneities  depend  not  only  on 
the  physical  properties  of  the  heterogeneity  (Wu  &  Aki,  1985),  but  also  need  include  the  vector 
quantity  k-a  to  reflect  the  fact  that  the  scattering  depends  not  simply  on  the  scalar  ratio  of  the 
’’size”  of  the  heterogeneities  to  the  wavelength  of  the  incident  wavefield,  but  also  on  the  spatial 
relation/orientation  of  the  incident  wavefieldfield  and  heterogeneities. 

In  future  modehng,  AN  ISO  heterogeneities  with  varying  aspect  ratios,  standard  deviations, 
physical  properties,  etc.,  will  be  incorporated  into  models  with  various  background  structures 
such  as  low  velocity  surface  layers,  and  surface  and  layer /half-space  boundary  topography.  Our 
modeling  has  also  highlighted  the  limitations  of  using  a  two-dimensional  routine  to  investigate  what 
is  truly  a  three-dimensional  phenomena.  Future  modeling  will  be  three-dimensional  and  include 
heterogeneities  with  ~  ay  >  a^. 
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